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ABSTRACT 

We  report  high  spatial  and  spectral  resolution  [Fe  ii]  1.644  /im  Fabry-Perot  imaging  observations  of 
M82.  We  present  extinction-corrected  [Fe  ii]  images  and  discuss  the  nature  of  compact  [Fe  n]  emission 
regions  revealed  by  these  new  data.  We  conclude  that  these  [Fe  n]  sources  trace  a  population  of  super¬ 
nova  remnants  in  M82  that  are  substantially  older  than  those  revealed  previously  on  6  cm  radiographs. 
In  addition,  we  find  that  M82  contains  a  distributed  [Fe  ii]  emission  component  that  is  extended  along 
the  southern  minor  axis  and  that  accounts  for  90%  of  the  galaxy’s  [Fe  ii]  luminosity.  We  interpret  this 
extended  emission  as  tracing  disk  material  entrained  in  a  super  wind  that  has  broken  out  of  the  galactic 
disk  to  the  south.  We  find  that  the  [Fe  n]/Bry  line  ratio  throughout  M82  correlates  with  the  age  of  the 
starburst  as  reflected  by  the  color  of  the  photospheric  emission  from  the  galaxy’s  stars.  This  correlation 
suggests  that  the  [Fe  ii]  emission  regions  in  M82  are  colocated  with  a  post-main-sequence  stellar  popu¬ 
lation.  The  engineering  details  of  our  Fabry-Perot  imaging  methodology  are  also  discussed. 

Subject  headings:  galaxies:  individual  (M82)  —  galaxies:  stellar  content  —  galaxies:  structure  — 
infrared:  galaxies  —  supernova  remnants  —  techniques:  interferometric 


1.  INTRODUCTION 

The  possibility  that  infrared  [Fe  ii]  emission  can  reveal 
the  supernova  activity  of  galaxies  has  stimulated  increasing 
observational  work  over  the  past  few  years.  Issues  sur¬ 
rounding  [Fe  ii]  as  a  tracer  of  supernova  remnants  in  gal¬ 
axies  via  grain  processing  were  discussed  by  Greenhouse  et 
al.  (1991),  hereafter  Paper  I.  Since  that  time,  other  interpre¬ 
tations  of  [Fe  ii]  emission  have  been  proposed  (Mouri, 
Kawara,  &  Taniguchi  1993),  progress  has  been  made  in 
understanding  relationships  between  infrared  [Fe  n]  lines 
and  other  emission  components  (e.g.,  Forbes  &  Ward  1993; 
Colina  1993;  Lumsden  &  Puxley  1995),  new  imaging  spec¬ 
trophotometry  of  [Fe  ii]  in  nearby  galaxies  has  been 
obtained  (Forbes  et  al.  1993;  Moorwood  &  Oliva  1994;  van 
der  Werf  et  al.  1993 ;  Blietz  et  al.  1994),  and  new  calculations 
of  Fe+  energy  levels  and  collisional  excitation  rate  coeffi¬ 
cients  have  been  produced  (Pradhan  &  Zang  1993;  Zang  & 
Pradhan  1995 ;  Nahar  1995). 

Emission-line  imaging  of  nearby  galaxies  has  revealed 
two  categories  of  near-infrared  [Fe  ii]  emission:  (1)  a  dis¬ 
tributed,  uniform  [Fe  n]  background  emission  and  also  (2) 
compact  [Fe  ii]  sources.  The  distributed  emission  has  been 
seen  to  correlate  spatially  with  global  nuclear  outflows  such 
as  starburst  super  winds  or  active  galactic  nucleus  radio 
jets.  The  compact  [Fe  ii]  sources  often  show  weak  or  no 
spatial  correlation  with  other  infrared  line  emission  such  as 
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H2,  H  i,  [Ne  ii],  or  continuum  emission  in  the  K  band  or  at 
6  cm. 

The  near-infrared  [Fe  ii]  image  data  obtained  to  date 
consist  of  low  spectral  resolution  narrowband  filter  images 
or  low  spatial  resolution  images  assembled  from  long-slit 
spectra.  In  this  paper,  we  present  the  first  near-infrared 
[Fe  ii]  Fabry-Perot  images  of  a  galaxy — a  technique  that 
yields  both  high  spectral  and  high  spatial  resolution.  In  the 
following  section,  we  discuss  our  methodology  for  Fabry- 
Perot  imaging,  the  image  analysis,  and  technical  details  of 
our  M82  [Fe  n]  observations.  In  §  3,  we  discuss  our  Fabry- 
Perot  images  of  M82  and  new  clues  that  they  reveal  about 
the  origin  of  M82’s  [Fe  ii]  emission. 

2.  OBSERVATIONS 

2.1.  Methodology :  Prefiltering  Facility  Cameras  with 
Fabry-Perot  Imaging  Systems 

Recent  advances  in  low  noise  detector  arrays  have  made 
Fabry-Perot  imaging  an  increasingly  practical  technique 
for  moderate-resolution  2/A2  ~  103  imaging  spectroscopy 
at  near-infrared  wavelengths.  High  system  accuracy  and 
reliability  can  be  achieved  with  commercially  available 
etalons  such  as  those  produced  by  Queensgate  Instruments 
LTD.8  Optimal  optical  coupling  between  the  Fabry-Perot 
and  camera  subsystems  is  critical  for  acquisition  of  data  on 
faint,  extended  objects  such  as  galaxies  and  has  not  been 
adequately  discussed  in  recent  articles  describing  infrared 
Fabry-Perot  imaging  systems.  As  a  consequence,  we  begin 
by  briefly  reviewing  engineering  details  that  govern  the 
optical  coupling  of  Fabry-Perot  etalons  to  infrared 
cameras.9  Details  of  our  M82  observations  and  data 
reduction  are  discussed  in  §  2.2. 

In  moderate  spectral  resolution  infrared  applications 

8  Queensgate  Instruments  LTD,  Queensgate  House,  Waterside  Park, 
Bracknell,  RG12  1RB,  England. 

9  For  a  more  complete  treatment  of  Fabry-Perot  formalism,  we  recom¬ 
mend  Vaughan  (1989)  and  Atherton  et  al.  (1981). 
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involving  relatively  high  f-ratio  (>f/15)  telescopes,  an 
optimal  optical  configuration  for  a  Fabry- Perot  camera 
system  is  one  in  which  the  etalons  are  located  directly  in  the 
telescope  beam.  The  standard  alternative  configuration 
incorporates  collimating  optics  that  add  complexity,  reflec¬ 
tion,  and  transmission  losses.  When  a  Fabry-Perot  is  used 
directly  in  a  converging  or  diverging  beam  of  cone  angle  cp, 
its  resolving  power  is  limited  to  Rheam  =  8/<p2.  The  total 
system  finesse  is  a  quadrature  sum  of  terms  due  to  etalon 
plate  reflectivity,  defects  of  the  plate  parallelism  and  surface 
figure,  as  well  as  geometrical  and  diffraction  optical  defects 
on  the  incident  beam.  Optimal  etalon  transmission  is 
achieved  when  the  reflectivity  finesse  is  much  less  than  each 
defect  finesse  term.  This  condition  can  often  be  met  at  slow 
Cassegrain  focal  ratios  without  employing  collimating 
optics.  For  example,  an  etalon  with  reflectivity  finesse  Nr  = 
50  operating  in  order  n  =  20  will  have  a  resolving  power  of 
R  =  1/A/,  ~  nNr  =  103.  A  f/27  telescope,  such  as  the  one 
used  here  (§  2.2),  yields  Rbeam  —  5.8  x  103  >  R  and  a  con¬ 
vergence  defect  finesse  /Vhcam  =  Rhciim/n  ~  290  >  Nr.  The 
corresponding  reduction  in  etalon  transmission  due  to 
beam  divergence  is  1  -  rbeam  =  1  -  (Nbeam//Vr)  arctan 
(Ar/Abeam)  =  0.01.  One  can  see  that  the  collimating  optics 
would  afford  no  significant  enhancement  of  this  etalon’s 
performance  and  would  increase  cost  and  complexity  while 
reducing  the  system  throughput  as  a  result  of  reflection  and 
transmission  losses.  The  former  loss  can  be  particularly  sig¬ 
nificant  when  the  collimating  optics  must  pass  a  broad 
wavelength  range  precluding  the  use  of  optimized  anti¬ 
reflection  coatings. 

Unnecessary  use  of  Fabry-Perot  collimating  optics  can 
result  in  nonmonochromatic  images  that,  in  infrared  appli¬ 
cations,  significantly  increase  observing  time.  For  example, 
a  Fabry-Perot  transmission  function  has  maxima  given  by 
10  =  2d  cos  9/n,  where  20  is  the  wavelength  of  the  nth  order 
maxima,  d  is  the  effective  gap  length,  and  9  is  the  incidence 
angle  of  the  chief  ray  from  a  point  of  interest  in  the  object 
field.  Chief  rays  arriving  from  field  points  across  the  field  of 
view  (FOV)  differ  in  incidence  angle  by  an  amount  that 
depends  on  the  optics  preceding  the  Fabry-Perot.  This 
effect  can  be  taken  into  account  in  image  processing  (§  2.2) 
if,  and  only  if,  an  image  data  cube  is  recorded  over  a  suffi¬ 
cient  number  of  gap  settings  to  fully  sample  an  emission  line 
over  the  entire  FOV. 

In  the  near-infrared,  the  background  emission  due  to  OH 
airglow  and/or  thermal  emission  is  generally  bright  and 
time  variable  on  the  scale  of  typical  narrowband  (R  =  103), 
background-limited  integration  times  ( ~  200  s).  As  a  result, 
explicit  sky  images  must  be  obtained  in  concert  with  the 
object  images  at  each  gap  setting.  A  sky-object-object-sky 
sequence  is  typically  used  to  remove  first-order  temporal 
variations  in  the  sky  emission.  Hence,  in  contrast  to  optical 
Fabry-Perot  applications,  each  extra  gap  setting  required 
by  the  above  9  dependence  results  in  at  least  four  extra 
integrations.  As  a  consequence,  in  infrared  applications,  it  is 
critically  necessary  to  minimize  the  range  of  chief  ray  inci¬ 
dence  angles  on  the  Fabry-Perot. 

The  range  of  incidence  angles  A 9  corresponding  to  a 
given  FOV  is  often  negligible  when  collimating  optics  are 
not  used.  In  this  case,  A 9  =  a/ j/,  where  a  is  the  field  stop 
diameter,  and  i jt  is  the  distance  between  the  focal  plane  and 
the  exit  pupil  of  the  telescope  (for  a  Cassegrain  telescope, 
\\i  =  primary  mirror  focal  length  +  back  focal  distance, 
e.g.,  see  Schroeder  1987).  The  focal  length  of  a  collimating 


lens  group  that  yields  similar  A 9  is  often  difficult  and  expen¬ 
sive  to  package,  particularly  in  cryogenic  applications. 

2.2.  M82  [Fe  n]  Observations  and  Data  Reduction 

In  Paper  I,  we  obtained  small  aperture  (277)  grating 
spectra  of  two  compact  nonthermal  continuum  sources  in 
M82.  We  speculated  that  high  spatial  resolution  Fabry- 
Perot  imaging  could  further  reveal  relationships  between 
[Fe  n]  emission  and  supernova  remnants  in  M82  and  other 
galaxies.  Here  we  report  on  Fabry-Perot  imaging  obser¬ 
vations  of  the  a4  D7/ 2  -» a4  Fg/2  (1.644  /an)  Fe+  forbidden 
line  in  M82  obtained  to  further  explore  the  origin  of  this 
emission. 

The  data  were  obtained  with  the  NASM/NRL  Fabry- 
Perot  system  and  the  University  of  Rochester  Third  Gener¬ 
ation  Camera  at  the  Wyoming  Infrared  Observatory 
(WIRO)  2.34  m  telescope  (f/27)  during  1991  May  12  UT.  A 
single  spatial  position  centered  on  the  nucleus  of  M82  was 
observed  using  a  2562  InSb  SBRC  focal  plane  array  with  a 
plate  scale  of  0'.'25  pixel-1.  The  spatial  resolution  was 
limited  by  seeing  to  1''3  FWHM  or  20  pc  at  the  distance  of 
M82  (3.2  Mpc;  Tammann  &  Sandage  1968).  All  of  the 
images  shown  here  have  been  block  averaged  to  a  display 
scale  of  0'.'5  pixel-1.  Image  rotation  was  determined  at  the 
telescope,  and  the  images  are  shown  here  with  pixel 
columns  oriented  north-south  to  within  approximately  O'.'l. 
Coordinate  axes  are  shown  in  units  of  arcseconds  relative  to 
the  K-band  nucleus:  a(1950)  =  9h5lm43:47,  5(1950)  = 
69°55'00'.'3  (Pipher  et  al.  1987). 

Four  Fabry-Perot  data  cubes  were  obtained  on  the 
[Fe  ii]  1.644  line  with  a  spectral  resolution  of  375  km  s-1. 
Each  data  cube  consists  of  a  three-dimensional  image  array 
(a,  d,  z),  where  a  and  d  are  coordinates  on  the  sky  and  z  is 
proportional  to  the  etalon  gap  length.  The  M82  data  cubes 
consist  of  nine  Fabry-Perot  gap  settings  chosen  to  sample 
the  line  profile  in  steps  of  (1/3)R  and  the  underlying  contin¬ 
uum  at  (5/3 )R  on  each  side  of  line  center.  At  each  gap 
setting,  four  images  were  recorded  following  a  standard  sky- 
object-object-sky  photometry  sequence,  thus  yielding  a 
total  data  set  of  144  images.  These  images  were  individually 
linearized,  flat-fielded,  sky-subtracted,  spatially  registered, 
and  calibrated  using  standard  infrared  image  reduction 
techniques.  Flux  calibration  was  determined  relative  to  the 
star  BS  3888  (F0  iv).  Spatial  registration  was  determined  by 
cross-correlation  analysis  on  the  point  source  K6  (Dietz  et 
al.  1986),  which  is  bright  in  [Fe  n]  and  off-line  continuum. 

With  the  above  reduction  complete,  the  remaining  task  is 
to  extract  continuum-subtracted  emission-line  images 
and/or  velocity  maps  from  the  processed  data  cubes.  Two 
optical  effects  must  be  taken  into  account  in  this  extraction : 
(1)  the  etalon  plates  typically  have  figure  defects,  and  (2)  the 
data  cubes  are  curved  because  of  the  finite  range  of  chief  ray 
incidence  angles  A 9.  Queensgate  Instruments  LTD  model 
ET50WF  etalons  represent  the  state  of  the  art  in  the  near- 
infrared  and  consist  of  water-free,  fused  silica  plates  with 
reflective  dielectric  coatings.  The  coatings  are  applied  at  an 
elevated  temperature,  and  a  large  number  of  layers  are 
required  to  achieve  broadband  wavelength  coverage.  The 
expansion  coefficient  mismatch  between  the  coating  and  the 
plate  substrate  produces  a  large  surface  stress  upon  cooling 
to  room  temperature,  resulting  in  a  bowl-shaped  surface 
deformation.  For  near-infrared  etalons  with  a  fractional 
wavelength  coverage  of  0.25  and  a  clear  aperture  of  50  mm, 
the  best  commercially  available  etalons  typically  exhibit  a 
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Camera  Lens 


Fig.  1. — Optical  layout  of  the  Fabry-Perot  camera  system  used  here 


1.642  1.644  1.646  1.648 


1.642  1.644  1.646  1.648 


flatness  defect  of  at  least  3  nm  edge  to  center  (2/200  at 
2  =  633  nm),  resulting  in  an  etalon  edge-to-center  gap 
variation  of  roughly  6  nm. 

For  etalons  used  in  a  collimated  beam,  this  gap  variation 
reduces  the  effective  finesse.  However,  in  the  optical  con¬ 
figuration  used  here  (Fig.  1),  the  radially  symmetric  defor¬ 
mation  produces  a  radial  dispersion  across  the  array  FOV. 


1.642  1.644  1.646  1.648 


1.642  1.644  1.646  1.648 


1.642  1.644  1.646  1.648  1.642  1.644  1.646  1.648 

Fig.  2. — Spectra  extracted  from  the  Fabry-Perot  data  cubes  described  in  §  2.2  at  the  spatial  center  of  [Fe  n]  sources  listed  in  Table  1.  A  typical  1  a  error 
bar  is  shown  in  the  upper  left-hand  corner  of  the  top  left  panel. 
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TABLE  1 

M82  [Fe  ii]  1.644  /an  Compact  Sources 


Offsets  from  Nucleus1* 


Fig.  4 a  Fig.  4 b  FWHMb  -  -  N( H) 

Source  Association  (west,  south)  (west,  south)  (pc)  Av°  Observed  Intrinsic1*  Observed  Intrinsic**  (ppm) 


Fe  1 .  ...  4.2, 4.2  3.7,3 .3  43  9  12.5  ±  0.50  49.60  24  ±  2  189  2.36 

Fe  2 .  ...  6.5, 0.9  6.5, 0.9  37  9  8.55  ±  0.04  33.50  17  ±  2  147  1.37 

Fe  3 .  ...  11.2,1.4  11.2,1.4  22  6  2.49  ±  0.03  6.59  7.3  ±  0.1  39.6  1.50 

Fe  4 .  ...  ...  10.2,4.2  37  9  ...  27.1  ...  150  1.63 

Fe  5 .  K56  23.7,7.3  23.7,7.3  22  ...  1.57  ±  0.05  ...  2.09  ±  0.4 

Fe  6 .  K6e  28.8,10.2  28.8,10.2  22  ...  1.54  +  0.06  ...  2.4  ±  0.4 


11  Offsets  in  arcseconds  relative  to  a(1950)  =  9h51m43147,  <5(1950)  =  69°55'00'.'3. 
b  FWHM  source  diameter. 

**  Extinction  data  from  Satyapal  et  al.  (1995)  averaged  over  listed  source  diameter. 
d  10“ 17  W  m“ 2  integrated  over  listed  source  diameter. 

**  Corrected  for  extinction. 

f  [JV(Fe)/JV(H)]0  =  47  ppm  (Grevesse  &  Anders  1989). 
g  Dietz  etal.  (1986). 


The  etalons  used  to  obtain  the  observations  reported  here 
have  a  plate  flatness  of  1/100,  yielding  a  wavelength  shift  of 
1  nm  between  the  center  and  the  edge  of  the  array  FOV.  As 
noted  in  §  2.1,  a  radial  dispersion  can  also  be  introduced  by 
the  telescope  optics.  In  the  optical  configuration  of  Figure  1, 
the  detector  array  is  the  effective  field  stop  a  =  11.3  mm, 
and  at  WIRO,  \j/  =  5.89  m.  Thus,  Af)  =  0.002  rad,  yielding  a 
negligible  wavelength  shift. 

At  the  telescope,  we  measure  the  actual  shift  due  to  the 
combination  of  these  effects  by  obtaining  data  cubes  on  a 
monochromatic  flat  field  produced  by  scattering  discharge 
lamp  light  off  of  a  diffusing  screen.  An  Airy  function  is  then 
fitted  to  the  lamp  line  at  each  of  2562  spatial  pixels  in  the 
lamp  data  cube.  The  spatial  variation  of  the  effective  etalon 
gap  is  thus  measured,  and  the  basis  for  a  relative  wave¬ 
length  calibration  of  each  spatial  pixel  in  the  astronomical 
data  cubes  is  obtained.  This  procedure  ensures  that  extrac¬ 
tion  of  astronomical  line  images  and  velocity  maps  are  not 
affected  by  the  above  optical  effects.  Absolute  wavelength 
calibration  is  obtained  by  observation  of  lamp  lines  closely 
bracketing  the  astronomical  line  of  interest  and  from  OH 
airglow  lines  in  the  sky  position  data  cubes. 

The  gap  space  data  cubes  (a,  d,  z)  are  thus  converted  to 
wavelength  space  (a,  d,  2).  A  linear  continuum  baseline  and 
line  profile  function  are  then  fit  to  all  of  the  data  at  each 
(a,  <$).  We  typically  use  an  Airy  function  convolved  with  a 
Gaussian  as  a  profile-fitting  function.  If  the  width  of  the 
best-fit  convolution  exceeds  the  instrument-response  func¬ 
tion,  a  Gaussian  deconvolution  is  performed  to  recover  the 
true  profile  function  of  the  emission  line.  However,  in  the 
case  of  M82,  the  [Fe  ii]  emission  was  spectrally  unresolved 
at  all  (a,  <5)  observed.  Typical  spectra  extracted  from  the 
reduced  data  cubes  on  [Fe  ii]  sources  listed  in  Table  1  are 
shown  in  Figure  2. 

At  each  pixel,  a  linear  continuum  baseline  is  subtracted, 
and  the  Gaussian  line  profile  or  the  instrument-response 
function  is  integrated  over  the  bandwidth  of  the  order¬ 
sorting  filter.  The  result  is  high  spatial  resolution  images  of 
integrated  line  intensity,  observed  line  wavelength,  line 
width,  and  underlying  continuum  flux.  Performing  these 
calculations  independently  at  each  (a,  <5)  ensures  that  all  of 
these  quantities  are  determined  accurately  over  the  whole 
FOV. 


3.  M82  [Fe  n]  emission 
3.1.  Morphology 

Figure  3  (Plate  1)  (also  Fig.  4a  [PI.  2])  is  an  image  of 
[Fe  ii]  1.644  pm  integrated  line  intensity  in  M82.  The 
prominent  A -band  (2.2  pm)  continuum  sources  mapped  by 
Dietz  et  al.  (1986)  and  the  radio  sources  mapped  by 
Kronberg,  Biermann,  &  Schwab  (1985)  are  plotted  on  the 
image  for  reference.  Two  general  categories  of  emission  are 
seen:  (1)  a  distributed  [Fe  ii]  background  emission  extend¬ 
ing  throughout  the  nuclear  region  of  the  galaxy  and  (2) 
compact  [Fe  ii]  sources  at  several  locations  southwest  of 
the  nucleus.  Compact  sources  were  identified  using  the  task 
DAOFIND  in  the  DAOPHOT  software  package  in  IRAF 
(Stetson  1987).  Six  compact  sources  were  identified  and  are 
labeled  on  Figure  4.  Their  positions,  FWHM  source  diam¬ 
eter,  and  integrated  line  flux  are  listed  in  Table  1. 

There  is  excellent  spatial  correspondence  between  the 
two  unresolved  [Fe  ii]  sources  in  the  southwest  corner  of 
the  image  and  the  Dietz  et  al.  sources  K5  and  K6.  These 
sources  are  strong  [Fe  ii]  point  sources  that  have  not  been 
detected  in  the  radio  continuum  (e.g.,  Huang  et  al.  1994; 
Muxlow  et  al.  1994;  Carlstrom  &  Kronberg  1991; 
Kronberg  et  al.  1985).  A  ridge  of  [Fe  n]  emission  is  associ¬ 
ated  with  K5  and  extends  to  K6.  The  eastern  terminus  of 
this  ridge  corresponds  to  a  point  of  high  foreground  extinc¬ 
tion  mapped  by  Satyapal  et  al.  (1995),  and  we  suspect  that 
this  ridge  actually  extends  eastward  along  the  disk  of  the 
galaxy,  connecting  with  the  [Fe  ii]  emission  seen  in  the 
inner  150  pc  region  of  the  nucleus  (Kl)  where  several 
compact,  spatially  resolved  [Fe  n]  sources  are  seen. 

In  Table  1,  we  list  observed  line  fluxes  and  positions  for 
the  prominent  [Fe  ii]  sources,  and  also  include  Pa/1  and  Av 
data  from  Satyapal  et  al.  (1995).  The  listed  fluxes  were  spa¬ 
tially  integrated  over  circular  regions  of  diameter  equal  to 
the  FWHM  extent  of  the  sources.  We  find  that  K5  and  K6 
(Fe  5  and  Fe  6  in  Table  1)  exhibit  a  [Fe  n]  1.644/Pa/i  ratio 
that  is  significantly  higher  than  any  other  [Fe  n]  source  in 
M82.  One  would  infer  a  [Fe  ii]  1.644/Bry  ratio  of  approx¬ 
imately  4,  or  a  [Fe  ii]  1.257/H/J  ratio  of  0.2  for  the  source 
K6.10  The  former  value,  and  the  upper  limit  to  the  source 
size  imposed  by  our  spatial  resolution,  suggest  that  K6  is 
comparable  in  size  and  emission  to  the  brightest  [Fe  ii] 
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TABLE  2 


M82  [Fe  n]  and  Pa/?  Photometry  (10  17  W  m  2)  Compared  with  Other  Data 


Source 

Beam 

(arcsec) 

A  a 

A-V 

This  Paper 

Paper  I 

Lester  et  al. 
([Fe  ii]  1.6) 

[Fe  ii]  1.6 

[Fe  ii]  1.3b 

Pa  /?c 

[Fe  n]  1.3 

Pa/? 

Nucleus . 

3.8 

5.0 

13.9 

16.0 

K2 . 

2.7 

6.0 

9.63 

4.60 

3.10 

4.60 

3.89 

41.9  +  58.1 . 

2.7 

6.0 

8.95 

4.80 

2.85 

13.5 

5.55 

44.0  +  59.5 . 

2.7 

3.2 

8.50 

6.50 

2.61 

23.0 

5.11 

3  Extinction  data  from  Satyapal  et  al.  (1995)  averaged  over  listed  aperture. 
b  Derived  from  [Fe  n]  1.644  using  listed  Av. 
c  Pa/?  data  from  Satyapal  et  al.  (1995). 


sources  in  NGC  253  (Forbes  et  al.  1993).  The  latter  value  is 
in  agreement  with  model  predictions  of  [Fe  n]  1.257/H/l 
behind  100  km  s” 1  grain-processing  shocks  (Shull  &  Draine 
1987). 

In  Table  2,  we  present  [Fe  n]  and  Vafl  fluxes  at  positions 
and  in  apertures  chosen  to  facilitate  comparison  with  M82 
[Fe  n]  data  published  by  others.  We  find  generally  good 
correspondence  with  results  reported  in  Paper  I  and  Lester 
et  al.  (1990)  on  the  nucleus  and  K2.  However,  we  find  lower 
[Fe  ii]  and  Pa/8  fluxes  than  reported  in  Paper  I  on  the  radio 
sources  44.0  +  59.5  and  41.9  +  58.1.  In  Paper  I,  we  assumed 
that  the  [Fe  ii]  surface  brightness  was  distributed  uniformly 
within  the  2"7  diameter  circular  aperture  of  the  grating 
spectrometer.  However,  these  new  high  spatial  resolution 
Fabry-Perot  images  reveal  that  this  emission  has  spatial 
structure  on  a  smaller  scale.  As  discussed  in  §  3.3,  we  find 
that  the  [Fe  ii]  fluxes  associated  with  the  radio  sources 
44.0  +  59.5  and  41.9  +  58.1  in  Paper  I  were  dominated  by 
compact  sources  on  the  periphery  of  the  beam. 

One  expects  the  [Fe  ii]  morphology  shown  in  Figure  3  to 
be  affected  by  nonuniform  foreground  extinction.  Satyapal 
et  al.  (1995)  derived  a  foreground-extinction  map  of  M82 
from  Fabry-Perot  images  of  Bry  and  Pa/f  obtained  with  the 
same  observing  system  at  the  same  spatial  and  spectral 
resolution  as  the  [Fe  ii]  data  presented  here.  We  have  used 
this  infrared-extinction  map  to  correct  Figure  4a,  assuming 
that  no  variable  extinction  exists  between  the  [Fen]  and 
the  H  recombination  line-emitting  regions.  Combination  of 
the  statistical  uncertainty  in  the  Av  and  [Fe  ii]  data 
increases  the  noise  in  the  resulting  corrected  image.  Those 
regions  where  the  extinction-corrected  [Fe  ii]  flux  was 
determined  with  >2  a  confidence  are  shown  in  Figure  4b 
(Plate  2).  Information  is  lost  on  the  K5  and  K6  region  of 
Figure  4b,  but  a  high  signal-to-noise  ratio  remains  through¬ 
out  the  inner  150  pc  region  that  contains  sources  Fe  1-4. 
The  source  Fe  4  is  revealed  with  high  contrast  by  this 
procedure. 

One  can  see  that  these  sources  are  intrinsically  compact 
and  do  not  result  from  holes  in  a  foreground-extinction 
screen.  The  brightest  source  is  located  just  west  of  the 
nucleus  and  consists  of  a  110  pc  x  40  pc  emission  region 
elongated  along  P.A.  =  40°,  with  maximum  emission  at  its 
southern  end  and  a  secondary  maximum  northwest  of  this 
position.  We  interpret  this  structure  as  two  sources  (Fe  1 
and  Fe  2,  Table  1)  that  exhibit  a  FWHM  diameter  of 

10  We  note  that  [Fe  n]  1.257/[Fe  ii]  1.644  =  1.36  (no  extinction), 
W/Pa/s  =  6.1,  and  =  5.7  (case  B,  104  K;  Hummer  &  Story  1987),  and 
we  caution  that  significant  departures  from  case  B  predictions  for  IBry/IHp 
are  often  seen  in  supernova  remnants  and  galactic  nuclei  (e.g.,  Oliva, 
Moorwood,  &  Danziger  1989;  Moorwood  1989). 


roughly  40  pc.  Source  Fe  4  (Table  1)  is  40  pc  x  30  pc 
FWHM  elongated  north-south.  We  note  that,  if  the  emis¬ 
sion  region  encompassing  Fe  1  and  Fe  2  were  interpreted  as 
resulting  from  two  adjacent  supernova  remnants,  and  if 
source  Fe  4  were  due  to  a  single  remnant,  these  remnants 
would  individually  be  similar  in  size  to  the  Cygnus  Loop 
(diameter  ~40  pc,  age  ~  5  x  104  yr).  In  contrast,  remnants 
corresponding  to  the  spatially  unresolved  sources  Fe  3,  Fe 
5,  and  Fe  6  would  have  to  be  at  least  a  factor  of  2  smaller 
and,  as  a  consequence,  much  younger  than  the  40  pc 
sources. 

We  note  that  no  strong  spatial  correlation  exists  between 
the  [Fe  n]  sources  and  the  [Ne  n],  [Ar  m],  [S  iv],  and  Bra 
sources  mapped  by  Achtermann  &  Lacy  (1995,  hereafter 
AL95).  The  sources  Fe  2  and  Fe  3  lie  roughly  15  pc  north  of 
the  AL95  sources  W1  and  W2.  We  find  no  [Fe  ii]  counter¬ 
part  to  the  AL95  source  El,  and  no  evidence  for  local 
maxima  in  the  [Ne  ii]  surface  brightness  at  the  position  of 
Fe  1,  Fe  4,  Fe  5,  or  Fe  6  is  apparent  in  the  AL95  line  maps. 
This  noncorrelation  is  not  surprising  due  to  the  factor  of 
3  difference  in  ionization  potential  between  photo- 
disassociation  region  (PDR)  species,  such  as  Fe+,  and  the 
H  ii  region  species,  such  as  [Ne  n],  although  the  latter  is 
also  modeled  to  trace  fast  shocks  (e.g.,  Shull  &  Draine 
1987).  We  also  note  that  no  [Fe  ii]  source  was  found  at  the 
position  of  the  [Ni  ii]  source  reported  by  Rank  et  al.  (1995). 

Figure  4c  (Plate  2)  is  a  composite  image  consisting  of 
[Fe  ii]  data  shown  in  Figure  4a  (red)  and  Bry  data  from 
Satyapal  et  al.  (1995)  (green).  Yellow  is  a  mixture  of  these 
components.  One  can  see  that  the  distributed  [Fe  ii]  emis¬ 
sion  is  extended  to  the  south  relative  to  Bry  in  the  general 
direction  of  the  galaxy’s  minor  axis.  This  extended  emission 
accounts  for  90%  of  M82’s  [Fe  ii]  luminosity  (Table  3).  A 
similar  morphology  has  been  seen  in  NGC  4945 
(Moorwood  &  Oliva  1994)  in  which  [Fe  n]  1.644  fim  is 
extended  symmetrically  along  the  minor  axis  relative  to  Bry. 
However,  in  M82,  the  [Fe  n]  emission  is  very  asymmetric, 
with  no  extension  to  the  north  relative  to  Bry,  and  no 
compact  [Fe  n]  sources  east  or  north  of  the  nucleus  are 
seen.  This  asymmetry  is  not  seen  in  the  1.64  gm  continuum. 


TABLE  3 

M82  [Fe  ii]  1.644  /im  Luminosities3 


Region 

Observed 

Extinction  Corrected 

30"  diameter  field . 

0.55 

1.20 

Total  of  compact  sources . 

0.03 

0.14 

Total  extended  emission . 

0.52 

1.06 

3  1040  ergs  s  1,d  =  3.2  Mpc. 
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Hence,  foreground  extinction  to  the  northern  half  of  the 
disk  is  not  a  possible  explanation.  We  note  that  this  north- 
south  asymmetry  is  also  seen  in  the  X-ray  continuum 
(Bregman,  Schulman,  &  Tomisaka  1995)  and  in  the  dis¬ 
tribution  of  6  cm  radio  continuum  point  sources  (Muxlow 
et  al.  1994). 

This  asymmetric  [Fe  ii]  spatial  distribution  most  likely 
traces  low  iron  depletion  gas  from  the  galactic  disk  that  has 
been  entrained  in  a  supernova-driven  wind.  The  north- 
south  asymmetry  suggests  that  the  wind  has  broken  out  of 
the  disk  in  a  southerly  direction  only.  One  could  also  specu¬ 
late  that  the  asymmetric  morphology  of  this  outflow  and 
the  southern  concentration  of  supernova  remnants  are 
related  to  the  H  i  stream  that  extends  south  from  M82  and 
bridges  the  M81/M82  system  (e.g.,  Gottesman  &  Welia- 
chew  1977;  Yun,  Ho,  &  Lo  1993). 

3.2.  Iron  Depletion  in  M82 

The  gas-phase  abundance  of  Fe+  can  be  derived  from  our 
measured  [Fe  ii]/Pa/7  ratios  as 


JV(Fe+) 
N(  H) 


9  x  10"6 


/([Fe  ii]  1.644) 

W) 


(1) 


(e.g.,  van  der  Werf  et  al.  1993),  where  we  have  assumed  a 
volume  emission  coefficient  for  Pa /?  of  jP^/n^  =  2.42  x 
10” 27  ergs  s-1  cm-3  sr-1  for  case  B  recombination  at 
T  =  104  K  (Table  4.2  in  Osterbrock  1974).  We  assume  low 
density  (ne  <  104  cm-3)  for  the  [Fe  n]-emitting  regions 
such  that  the  a4  D  term  of  Fe+  is  not  thermalized  (Paper  I), 
and  we  take  /V(Fe+)/A,r(Fe)  as  order  unity  in  shock  recombi¬ 
nation  zones  (Paper  I). 

From  equation  (1)  and  Table  1,  we  find  a  gas-phase  iron 
abundance  /V(Fe)//V(Fe)G  >  0.04  averaged  over  sources  Fe 
1-4,  and  on  sources  Fe  5  and  Fe  6,  we  find  A(Fe)/JV(Fe)0  > 
0.12.  For  M82  as  a  whole  (30"  diameter  field),  we  find 
L[Feiii  =  5.5  x  1039  ergs  s_1  (Table  3)  and  LBty  =  6.5 
x  1039  ergs  s_1  (Satyapal  et  al.  1995)  yielding  a  global 
gas-phase  abundance  of  A(Fe)/lV(Fe)0  >0.16.  The  lower 
limits  derive  from  our  assumption  that  all  of  the  iron  is 
singly  ionized.  These  data  show  that  the  iron  depletion  in 
these  sources  is  at  least  an  order  of  magnitude  lower  than 
typical  interstellar  values  along  reddened  lines  of  sight  in 
the  Galactic  interstellar  medium  (de  Boer,  Jura,  &  Shull 
1987),  and  is  at  most  roughly  equivalent  to  model  predic¬ 
tions  of  iron  depletion  behind  moderate-density  100  km  s_  1 
shocks  (McKee,  Chernoff,  &  Hollenbach  1984),  consistent 
with  the  idea  that  the  [Fe  ii]  morphology  in  M82  traces 
regions  of  grain  destruction  (Paper  I).  We  note  that  the 
global  L[Fell]/LBry  ratio  we  derive  is  consistent  with  radio 
hypernovae  models  of  Colina  (1993)  that  involve  a  Salpeter 
(a  =  2.35)  initial  mass  function  slope  with  an  upper  mass 
limit  of  25  M0 .  However,  L[Fe  n]  and  LBry  are  individually  a 
factor  of  2  in  excess  of  this  prediction.  The  compact  mor¬ 
phology  of  M82’s  [Fe  ii]  sources  strongly  suggests  that 
these  sources  are  supernova  remnants.  However,  at  the  high 
spatial  resolution  of  these  new  data,  we  find  that  these 
sources  are  distinct  from  the  population  of  radio  point 
sources  revealed  on  6  cm  radiographs  of  M82. 

3.3.  [Fe  ii]  and  Radio  Continuum  Sources 

In  Paper  I,  we  obtained  [Fe  ii]  1.257  spectra,  with  a 
spatial  resolution  of  42  pc,  of  two  nonthermal  6  cm  contin¬ 
uum  sources,  41.9  +  58.1  and  44.0  +  59.5  (Kronberg  et  al. 

1985) ,  and  a  K-band  continuum  source  K2  (Dietz  et  al. 

1986) .  We  compared  the  [Fe  ii]  emission  of  the  radio 


sources  with  K2,  and  also  with  archival  [Fe  ii]  data  on  K1 
reported  by  Lester  et  al.  (1990).  We  concluded  that  the 
[Fe  ii]  surface  brightness  had  local  maxima  at  the  position 
of  these  two  radio  sources.  At  the  higher  spatial  resolution 
(20  pc)  of  the  new  data  reported  here  (Fig.  4),  we  find  that 
the  [Fe  ii]  source  previously  attributed  to  44.0  +  59.5  is 
displaced  one  spatial  resolution  element  south  of  this  radio 
source.  We  also  find  that  the  [Fe  n]  emission  attributed  to 
41.9  +  58.1  is  due  to  emission  on  the  northeast  and 
southwest  periphery  of  the  2" 7  beam  used  in  our  past  mea¬ 
surements  (Paper  I),  such  that  the  [Fe  ii]  emission  at  this 
position  is  actually  attributed  to  sources  Fe  2  and  Fe  4 
rather  than  to  the  radio  source. 

Although  several  of  the  [Fe  ii]  sources  in  Table  1  are  in 
close  proximity  to  one  or  more  radio  continuum  sources, 
we  find  no  strong  spatial  correlation  between  the  [Fe  ii] 
and  the  radio  continuum  sources  mapped  by  Kronberg  et 
al.  (1985)  or  more  recently  by  Muxlow  et  al.  (1994).  In  terms 
of  total  spatially  integrated  luminosity,  [Fe  ii]  and  6  cm 
emission  components  are  tightly  correlated  (e.g.,  Forbes  & 
Ward  1994).  However,  in  terms  of  spatial  morphology, 
these  components  are  not  correlated.  Since  the  [Fe  ii]  point 
sources  make  only  a  small  direct  contribution  to  the  total 
[Fe  ii]  luminosity  at  these  wavelengths  (Table  3),  we  con¬ 
clude  that  this  luminosity  correlation  is  tracing  the  distrib¬ 
uted  components.  The  lack  of  spatial  correlation  between 
the  [Fe  n]  and  radio  point  sources  and  the  size  of  the 
spatially  resolved  [Fe  n]  sources  suggests  that  these  com¬ 
ponents  are  tracing  separate  populations  of  supernova  rem¬ 
nants  that  differ  significantly  in  age. 

The  radio  point  sources  in  M82  are  interpreted  as  young 
supernova  remnants  rather  than  radio  supernovae,  based 
on  the  temporal  development  of  their  6  cm  flux.  Many  of 
these  remnants  have  been  spatially  resolved.  Their  age  is 
constrained  by  their  size.  Muxlow  et  al.  (1994)  report  2.4  pc 
as  a  mean  remnant  diameter  for  the  population  and  derive  a 
corresponding  mean  age  of  ~  200  yr  assuming  preadiabatic 
expansion  at  5000  km  s_1.  In  contrast,  we  find  that  the 
[Fe  ii]  sources  are  typically  a  factor  of  20  larger  than  the  6 
cm  sources.  The  proportionality  between  remnant  diameter 
and  age  r  oc  t"  changes  as  the  remnant  evolves  through  free 
expansion  (n  =  1),  adiabatic  (n  =  §),  and  radiative  (n  =  5) 
phases  of  expansion  (e.g.,  Spitzer  1978).  We  note  that  the 
size  contrast  between  the  [Fe  11]  and  the  6  cm  remnants  in 
M82  is  similar  to  that  exhibited  by  the  Cygnus  Loop  (age 
5000  yr)  and  Cas  A  (age  300  yr).  A  natural  interpretation  for 
the  compact  [Fe  11]  sources  in  M82  is  that  they  trace  a 
population  of  supernova  remnants  that  are  at  least  an  order 
of  magnitude  older  than  those  observed  in  the  radio  contin¬ 
uum.  We  note  that  the  size  distribution  of  the  [Fe  11] 
sources  spans  at  least  a  factor  of  2.  Assuming  that  these 
remnants  are  in  radiative  expansion,  they  must  span  a 
factor  of  10-20  in  age.  Hence,  these  data  suggest  that  [Fe  11] 
remains  an  important  coolant  in  supernova  remnants  for  at 
least  103  yr  after  outburst.  We  note  that  [Fe  11]  emission 
associated  with  the  2  pc  radio  sources  would  be  diluted  by 
our  20  pc  seeing-limited  resolution.  As  a  consequence,  the 
[Fe  11]  background  emission  observed  throughout  the 
nuclear  region  could  possibly  inhibit  their  detection.  Simi¬ 
larly,  6  cm  emission  associated  with  the  bright  20-40  pc 
[Fe  11]  sources  would  be  difficult  to  detect  against  the  6  cm 
continuum  background  (T.  W.  B.  Muxlow  1996,  private 
communication).  However,  we  note  that  the  sources  Fe  5 
and  Fe  6  would  be  good  candidates  to  search  for  a  radio 
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continuum  counterpart  since  they  are  the  most  compact 
[Fe  ii]  sources  and  occur  at  a  location  where  the  6  cm 
continuum  background  is  a  minimum. 

We  find  that  the  [Fe  ii]  luminosity  of  the  sources  listed  in 
Table  1  exceeds  that  expected  in  Galactic  or  LMC  super¬ 
nova  remnants  (1033-1036  ergs  s_1;  Oliva  et  al.  1989)  by 
factors  ranging  from  1  to  several  orders  of  magnitude.  A 
similar  affect  was  noted  by  Kronberg  et  al.  (1985),  who 
reported  6  cm  compact  source  luminosities  that  exceed  Cas 
A  by  2  orders  of  magnitude.  Averaging  over  observed  data 
listed  in  Table  1,  we  estimate  the  typical  [Fe  n]  luminosity 
of  a  M82  supernova  remnant  as  at  least  1038  ergs  s_1  (cf. 
Colina  1993).  Assuming  that  this  emission  remains  roughly 
constant  as  the  remnant  ages,  we  find  that  the  total  [Fe  ii] 
luminosity  of  M82  can  be  accounted  for  by  roughly  102 
supernova  remnants.  If  their  [Fe  n]  emission  persists  for 
103  yr,  then  these  data  would  suggest  a  supernova  rate  of  at 
most  0.1  yr” 1,  consistent  with  that  derived  by  Van  Buren  & 
Greenhouse  (1994)  from  radio  data  on  M82. 

3.4.  [Fe  ii]  and  Photospheric  Continuum  Emission 

If  the  [Fe  n]/Bry  flux  ratio  is  tracing  the  relative  contri¬ 
bution  of  supernova  remnants  and  Lyman  continuum- 
producing  stars  within  a  stellar  cluster,  then  these  data  can 
yield  a  general  age  indicator  of  the  central  starburst  in  M82 
(Paper  I).  At  the  onset  of  a  starburst,  the  stellar  population 
should  be  characterized  by  a  large  H  recombination  line 
luminosity,  a  small  [Fe  ii]  luminosity,  and  blue  near- 
infrared  continuum  colors.  As  the  most  massive  stars 
formed  during  the  burst  evolve  off  the  main  sequence,  the 
[Fe  ii]  luminosity  will  increase  with  the  supernova  rate, 
emission  from  low-depletion  PDRs  will  increase  relative  to 
H  ii  regions,  and  the  photospheric  continuum  will  become 
redder. 

We  find  a  correlation  between  the  J—H  color  and  the 
extinction-corrected  [Fe  n]/Bry  flux  ratio  within  the  central 
500  pc  of  M82  (Fig.  5).  The  J  and  H  broadband  data  shown 
in  Figure  5  have  been  corrected  for  extinction  and  free-free 
and  free-bound  continuum  emission.  We  find  that  dust 
emission  and  our  measured  Pa/1  and  [Fe  n]  line  fluxes  have 
a  negligible  effect  on  the  J  —  H  color  of  M82  (see  Satyapal  et 
al.  1995,  1996).  As  a  result,  we  assume  that  the  corrected 
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Fig.  5. — The  J  —  H  color  of  photospheric  continuum  emission  vs. 
[Fe  ii]  1.644/Bry  line  ratio  within  the  central  500  pc  of  M82. 


color  shown  in  Figure  5  reliably  traces  direct  photospheric 
emission  from  M82’s  stars  (cf.  Koornneef  1983).  We  find 
that  the  [Fe  n]/Bry  line  ratio  throughout  M82  correlates 
with  the  age  of  the  starburst  as  reflected  by  the  color  of 
photospheric  emission  from  the  galaxy’s  stars.  This  corre¬ 
lation  suggests  that  the  [Fe  n]  emission  regions  in  M82  are 
colocated  with  the  post-main-sequence  stellar  population, 
and  lends  credence  to  the  idea  that  infrared  [Fe  n]/H  ii  line 
ratios  can  yield  a  relative  age-dating  technique  for  starburst 
clusters. 

3.5.  The  [Fe  n]  Rotation  Curve  of  M82 

Although  the  resolving  power  of  our  Fabry-Perot  etalons 
is  ^ 375  km  s-1,  our  data  reduction  algorithms  (§  2.2) 
enable  systemic  velocity  measurements  to  be  made  with  a 
precision  of  approximately  40  km  s_1  across  the  whole 
FOV.  As  a  consequence,  we  easily  resolve  the  rotation  of 
M82.  The  relative  [Fe  n]  emission-line  velocity  along  the 
major  axis  of  the  galaxy  is  shown  in  Figure  6.  Good  agree¬ 
ment  is  found  between  our  [Fe  n]  rotation  curve  and  those 
derived  from  Ha  measurements  (Gotz  et  al.  1990).  In  addi¬ 
tion,  our  velocity  gradient  is  in  close  agreement  with  [S  m] 
9069  A  measurements  by  McKeith  et  al.  (1993).  This  agree¬ 
ment  provides  confirmation  that  our  Fabry-Perot  obser¬ 
vation  and  data  reduction  techniques  (§  2.2)  are  accurate 
over  the  whole  FOV  shown  in  Figure  4. 

3.6.  Conclusions 

We  conclude  that  infrared  Fabry-Perot  imaging  is  a 
powerful  technique  for  high  spatial  resolution  imaging  spec¬ 
troscopy  in  applications  requiring  high  spectral  purity  rela¬ 
tive  to  circular  variable  filter  imaging,  and  wide  fields 
relative  to  those  obtained  by  integral  field  slit  spectro¬ 
meters.  Our  [Fe  ii]  Fabry-Perot  images  of  M82  have 
revealed  a  population  of  supernova  remnants  that  are  sub¬ 
stantially  older  than  those  revealed  on  6  cm  radiographs  of 
similar  spatial  resolution.  We  note  that  this  interpretation, 
combined  with  recent  evidence  for  radially  propagating 
starburst  activity  in  M82  (e.g.,  Satyapal  et  al.  1996),  would 
imply  that  the  [Fe  ii]  sources  should  be  more  centrally 
concentrated  than  the  6  cm  sources.  This  expectation  is 
consistent  with  observation. 
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Fig.  6. — The  [Fe  n]  1.644  /an  rotation  curve  along  the  major  axis  of 
M82.  The  distance  is  shown  relative  to  the  nucleus,  and  the  velocity  was 
derived  from  the  line  wavelength  measured  relative  to  the  nucleus  along 
P.A.  120°. 
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We  also  find  that  the  nuclear  region  of  M82  contains  a 
distributed  [Fe  n]  background  emission  that  accounts  for 
90%  the  galaxy’s  [Fe  n]  luminosity  and  traces  a  region 
where  the  gas-phase  abundance  of  iron  is  at  least  16%  of 
solar.  We  find  that  the  gas-phase  iron  abundance  correlates 
spatially  with  the  color  temperature  of  the  photospheric 
continuum,  suggesting  that  the  low  iron  depletion  regions 
are  associated  with  evolved  stars  rather  than  H  ii  regions. 

Among  the  most  interesting  aspects  of  M82’s  [Fe  ii]  mor¬ 
phology  is  its  asymmetry.  We  find  that  the  minor-axis 
extension  of  the  [Fe  ii]  emission  follows  the  north-south 
asymmetry  revealed  in  ROSAT  HRI  continuum  maps,  and 
no  compact  [Fe  ii]  sources  were  detected  east  or  north  of 
the  nucleus  where  a  [Ni  n]  source  was  recently  observed 
from  the  Kuiper  Airborne  Observatory.  Indeed,  at  20  pc 
resolution,  we  find  no  [Fe  ii]  compact  sources  in  spatial 
coincidence  with  young  (2  pc  diameter)  supernova  rem¬ 
nants.  Similarly,  no  6  cm  sources  of  sizes  comparable  to  the 
[Fe  ii]  sources  (20-40  pc  diameter)  have  been  detected. 

3.7.  Starburst  Fossil  Remnants 
Can  poststarburst  galaxies  be  distinguished  from  normal 
galaxies?  If  the  distributed  [Fe  n]  emission  observed  along 


the  minor  axis  of  M82  and  NGC  4945  is  tracing  low- 
depletion  disk  material  entrained  in  a  starburst  super  wind, 
then  after  the  cessation  of  starburst  activity,  this  material 
should,  in  many  cases,  occupy  a  halo  surrounding  the 
galaxy.  If  the  low-depletion  state  of  this  gas  were  persistent 
in  the  halo  environment,  then  it  could  provide  an  observ¬ 
able  fossil  remnant  of  the  starburst.  One  could  perhaps 
image  poststarburst  low-depletion  halos  via  low-excitation, 
potential  ground-state,  fine-structure  transitions  of  refrac¬ 
tory  ions  using  the  Stratospheric  Observatory  For  Infrared 
Astronomy  (SOFIA).  Examples  would  include  the  a6D( 7/2— 
9/2)  transition  of  [Fe  ii]  near  26  /on  and  the  a5D(3-4)  tran¬ 
sition  of  [Fe  hi]  near  23  /on.  Both  transitions  have 
low-excitation  potential  (550  and  628  K,  respectively)  and 
would  be  excited  in  a  quiescent  poststarburst  environment. 
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ERRATUM 

In  the  paper  “  Infrared  Fabry-Perot  Imaging  of  M82  [Fe  n]  Emission.  II.  Tracing  Extragalactic  Supernova  Remnants  ”  by 
Matthew  A.  Greenhouse,  Shobita  Satyapal,  Charles  E.  Woodward,  J.  Fischer,  K.  L.  Thompson,  W.  J.  Forrest,  J.  L.  Pipher, 
N.  Raines,  H.  A.  Smith,  D.  M.  Watson,  and  R.  J.  Rudy  (ApJ,  476,  105  [1997]),  there  are  errors  in  Figure  4  (Plate  2)  and  Table 
1.  In  Figure  4,  the  panels  of  the  plate  were  mislabeled  by  the  printer.  As  the  plate  appears,  the  top  panel  is  actually  Figure  4c 
and  the  bottom  panel  is  actually  Figure  4a.  In  Table  1,  diameters  given  in  the  fifth  column,  for  the  unresolved  sources  Fe  3, 
Fe  5,  and  Fe  6  should  appear  as  “  <  22  ”  rather  than  “  22  ”. 
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